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Dynamical problems of the spin-stabilized Telstar satellite, character- 
ized by spin decay, spin-precession damping, and spin-axis drift, are ana- 
lyzed in this paper. Both the eddy-current torques and the magnetic torques, 
which cause the above three phenomena, are evaluated. By extrapolation from 
the observed data, the characteristic time of the nearly exponential spin decay 
of the satellite is estimated to be about 330 days. A linear analysis of the 
precession damper is made, and the results are compared with experiments, 
showing that the satellite precession angle will diminish by a factor of e in a 
maximum time of 30 minutes. A qualitative description is given to illustrate 
the fundamental mechanism of spin-axis drift. Results of these analyses can 
be applied to any spin-stabilized satellite. 

I. INTRODUCTION 

For spin stabilization of a communications satellite, it is required that 
the satellite be statically and dynamically balanced so as to make the 
principal axis of maximum moment of inertia coincide with the axis of 
symmetry of the antenna pattern, about which the satellite is given an 
initial spin. This principal axis, referred to henceforth as the spin axis, is 
in line with the invariant angular momentum vector and is thus fixed in 
direction, as desired, in an inertial space, provided there are no external 
torques acting on the spinning satellite. However, as the satellite is spin- 
ning and traveling in the geomagnetic field, eddy-current and magnetic 
torques continuously act on the satellite so that the angular momentum 
changes its magnitude and direction, as characterized by spin decay and 
spin precession. As a consequence of spin decay, the satellite becomes less 
stable for the same external disturbing torques, and a tumbling motion 
may eventually result. The precession of the spin axis about the instan- 
taneous angular momentum vector will cause wobbling of the antenna 

2169 



2170 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1963 

pattern, indicating that the precession should necessarily be dissipated 
by means of a damping mechanism. Because of the continuous action of 
the torques, the angular momentum continuously changes its direction 
in the inertial space; meanwhile, the spin axis precesses about it and, due 
to the precession damping, tends to align with it. Thus, there results a 
gradual drift in direction of the spin axis (sometimes called long-term pre- 
cession), as already observed on the Telstar satellite. 

The above dynamics problems of the satellite — namely, spin decay, 
spin-precession damping, and spin-axis drift — are studied in this paper. 
In the discussion of spin decay, we will indicate the nature of the retard- 
ing torques resulting from eddy currents and magnetic hysteresis losses, 
analyze the observed spin decay phenomenon, and compute the 1/e char- 
acteristic time of the exponential decay. For the spin precession, a linear 
analysis of the precession damping mechanism will be given, and an ex- 
perimental comparison of the damping time will be outlined. Only a short 
descriptive analysis is given to the problem of spin-axis drift, since an 
exact evaluation of the rate and pattern of drift deserves a separate com- 
puter study. 

It is shown in the following that whenever a spinning rigid body under- 
goes energy losses (e.g., from internal friction) the axis of maximum 
moment of inertia or the spin axis, z, and the angular velocity, w, will 
tend to align with the angular momentum, J = * • <o, where 

<S> = IJx + I v yy + hzz 

(x, y, and z are the unit vectors along the principal axes) is the moment 
of inertia dyadic. If I z of the spin axis is the minimum, the spinning mo- 
tion is still stable; however, any energy dissipation will not reduce pre- 
cession but make the spin axis deviate away from J. A simple proof of 
the above is given in the following. The kinetic energy, E, of a spinning 
satellite with precession is written as 

2E = w«I>-g>. (1) 

By substituting J = 4» o> and <o = J<& _1 into the above, E can be ex- 
pressed in terms of the angular momentum, 

J = J(cos £6 + cos ijy + cos 6z), 

where cos £, cos r;, and cos are the direction cosines of J in the body co- 
ordinates 



2E 



= J-* -1 - J = J 2 (t cos 2 £ + i cos 2 v + y cos 2 d) (2) 
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or, as cos 2 = (1 — cos 2 £ — cos 2 77) 



2 *= j lMM) c ^ + fe4) cos% 



(3) 



If there exists energy dissipation at a rate assumed to be so slow as to 
produce no torque on the satellite, the torque-free rigid body motion of 
the satellite will tend toward the state of minimum energy. It is observed 
from (2) and (3) that the minimum kinetic energy state occurs at 
$ = 0° (or £ = V = 90°) for I z > I x , I y , or at = 90° (either £ ■ 
or 77 =0) for J z < L , h • This proves that in order to reduce the pre- 
cession angle, 0, by means of energy dissipation, l z of the spin axis should 
be the maximum. 



11. spin DECAY 

Spin decay results from energy losses in the form of both eddy currents 
and magnetic hysteresis when the satellite is spinning in the geomagnetic 
field. It is shown in the following that the hysteresis losses in the magnetic 
materials are much smaller than the eddy-current losses in the conducting 
materials at high spin rates simply because the geomagnetic field in the 
Telstar satellite orbit is relatively weak, ranging from 0.04 to 0.4 oersted.* 
Magnetic materials are contained in the nickel-cadmium cells of the 
battery, the magnetic shielding on circuit components, etc. No measure- 
ment of the magnetic hysteresis loops has been made on the components 
actually used in the satellite. However, measurements on similar com- 
ponents contemplated for use on a proposed satellite prior to the Telstar 
satellite have been made, based on which it was estimated (for different 
orbit parameters and a different spin-axis attitude from those of the 
Telstar satellite) that the time-average hysteresis loss is W = 1.6 ergs 
per cycle of rotation. It is believed that the above value can be used 
for a conservative estimate of spin-decay rate for the Telstar satellite 
because it contains less magnetic materials than had been anticipated. 
According to this value of W, the spin decay of the satellite due to hys- 
teresis losses alone is only about 1.5 rpm per year. 

Eddy currents are generated essentially in the following parts : (a) the 
aluminum shell of the electronics chassis, (b) the frames of square magne- 
sium tubing and equatorial antennas, and (c) the magnesium chassis 
frame assembly. An estimate of the eddy-current torque can be made if 
the electronics chassis is approximated as a thin spherical shell and the 



* These figures are based on a spherical harmonic representation of the geo- 
magnetic field with Hensen and Cain coefficients for the Epoch, 1960. 
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last two items in the above are approximated as circular loops of wire. 
The eddy-current torque acting on a thin spherical shell spinning at an 
angular velocity w can be shown to be 

Ti = Pl B X (B X ») (4) 

where B is the geomagnetic induction and pi = (2ir/3)a a d, with a = 
radius, d = thickness, and a = volume conductivity. The above expres- 
sion is correct only when the square of the nondimensional quantity, 
\naoud (n = free-space permeability), is negligibly small compared to 
unity, which is found to be true in the present case. The above torque 
can be resolved into two components, i.e., the component parallel to <o 

T\ } = -piBja, (5)* 

which tends to retard to(.Bj. = component of B normal to to) and the 
component normal to w 

T x = pifiiifixw (6) 

which contributes to the precession of the satellite (B\\ = component of 
B parallel to <■>). In the case of a circular loop, it can be easily shown 
that the time-average eddy-current torque, acting on a circular loop 
spinning about a diameter, tends only to retard g>; i.e. 

T 2 = —p 2 B ± 2 (o (7) 

where P2 = A J2R' (A = loop area, R' = l'/aA w = total resistance of 
the loop of wire, A w = cross-sectional area of the wire, I' = length of the 
loop). The above expression is correct only when the square of the non- 
dimensional quantity, coL'/R' (L' = inductance), is negligibly small 
compared to unity, which is found to be the case here. The other com- 
ponent normal to to has a zero time-average value. 

In general, the eddy-current retarding torque acting on a conducting 
body spinning in a magnetic field is proportional to B x 2 and w, or can be 
written approximately as 

T r = -pBju (8) 

where p is a constant and is determined by the conducting material and 
its geometry. Thus, for the Telstar satellite, if the electronics chassis is 
approximated as a thin spherical shell (of 9.5-inch radius and 0.1-inch 
thickness) and the frames are approximated as circular loops of wire, 
then the retarding torque, T T , acting on the satellite is the sum of T\ in 



* This expression is the same as that given in Ref. 1, p. 417, problem 12, after 
the square of the nondimensional quantity, f/xewwd, is neglected. 
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(5) and T 2 in (7) or p = pi + p 2 . It is calculated that pi = (584, 
p 2 = 25G, or p = 940 meter 4 /ohm in mks units. This value of p is of 
course too low, because many small conducting parts have not been con- 
sidered in the calculation. From the expressions of p\ and p 2 , one notices 
that pi is proportional to the fourth power of the radius of a spherical 
shell and p 2 to the square of the loop area. For this reason, the Telstar 
satellite was insulated at the equatorial antennas in such a way that the 
outer shell does not constitute a large continuous surface and that the 
frames do not form large continuous loops. 

The magnitude of p in (8) for the Telstar satellite can be measured by 
rotating a magnetic field normal to the spin axis while the angular deflec- 
tion of a torsion wire, which suspends the satellite along the spin axis, 
is recorded to determine the drag torque. Such an experiment has been 
devised by M. S. Glass and D. P. Brady. Measurements made on the 
prototype give p = 1355 meter 4 /ohm ±15 per cent. This measured 
value is believed to be somewhat high, because the magnetic field applied 
in the measurements was as high as 25 to 100 oersteds (at 23.4 rpm) in 
order to give significant angular deflection readings of the suspension 
wire; thus, the measured drag torque unavoidably includes losses due to 
full hysteresis loops described in the magnetic materials. In the actual 
case, the magnetic field along the orbit is only 0.04-0.4 oersted, and the 
losses due to minor hysteresis loops are much smaller. Besides, since the 
electromagnetic characteristics of the satellite may be different from one 
model to another, the value of p measured on the prototype may not be 
applied to the Telstar satellite with good accuracy. Nevertheless, it is 
believed that the value of p in meter 4 /ohm is bounded below by 940 and 
above by 1560. A later calculation by extrapolation from the observed 
data on the satellite showed that p is approximately equal to 1110. A 
further refinement of the evaluation of p might have been obtained from 
the instantaneous spin-decay rate which can be determined from the 
telemetry solar aspect data. However, as the obtained instantaneous 
spin-decay rate was too low to give any significant reading, such an at- 
tempt failed to yield any results. 

Because of proper functioning of the precession damper, the Telstar 
satellite is now spinning nearly about its principal z-axis of maximum 
moment of inertia. In this case, B± in (8) can be approximated as the 
component of B normal to the z-axis. Obviously, B ± is a function of time 
due to (a) the rotation of the slightly inclined geomagnetic field about 
the earth's spin axis, (b) the anomalies of the geomagnetic field, (c) 
the gradual drift in direction of the satellite spin axis, and (d) the varia- 
tion of orbital parameters due to the oblateness of the earth, notably the 
apsidal advance and the nodal regression (see Fig. 1). For the same rea- 
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sons, the magnetic hysteresis loss per cycle of rotation, W, is also a func- 
tion of time. The decay of spin rate due to both eddy-current and hys- 
teresis losses can be determined from the following equation 

I«> = -pBj(t)<a-l-W(t) (9) 

or, upon integration, 

« = exp(-f (p/I z )B ± 2 (t)dt) 

i (10) 

•[coo - £ 2^ exp (£ (p/I.)Bj(t) dty dt . 

where u> = u(t = 0). Let us now define day-average values BJ and 
Was 

BJ=\ f B x \t)dt (11) 



and 



f W(t) exp (7' (p/I,)BS(t) dt) dt 

W=- j f-5 '—. (12) 

£ mp(£ (p/I.)Bj{t)dt\dt 



Then (10) can be written as 

r , i w(t) 

CO = \ b)Q •+- =- 



e-" T 



Wit) 



2*pBAt)- ^pBJ(t) 



(13) 



where 



r = IJpBAD- (i4: 



The time £ in (13) is in units of days, and BJ(t) and W(t) are functions 
of t. Note that if the term W/2irpBJ, which is much smaller than w 
in the case of the Telstar satellite, is neglected in (13), the spin decay 
is exponential with time with, however, a time-dependent t. 

A plot of BJ{t) is given in Fig. 2 from the day of launch (July 10, 
1962) up to December 31, 1962. The day-average BJ(t) is obtained by 
taking the arithmetic mean of the time-average values of B ± ~ per pass 
for approximately nine passes a day. The latter values are computed* 



Computations were provided by J. D. Gabbe. 
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Fig. 2 — Spin rate and B ± 2 vs time. 



from a spherical harmonic representation of the geomagnetic field, taking 
into account the continuous variations of the spin-axis direction and of 
the orbital parameters. Because of these combined effects, the variation 
of B ± 2 is sinusoidal with time with, however, variable amplitude and 
period. The major contribution to this variation is believed to be due 
to the apsidal advance in the orbital plane. In order to see this, let us 
plot in Fig. 3 the time-average values of B±- per pass versus the geo- 
graphical longitude of the perigee in the beginning of the pass on the 
days of July 15, 18 and 21, 1962. The latitudes of the perigee on those 
three days were ±5° within the geographical equator. It is shown on 
these curves that Bx 2 is relatively low when the perigee falls in the region 
over South America where the geomagnetic field strength is depressed. 
The center of this region falls at approximately 25°S latitude and 45°W 
longitude. (See Ref. 2 for details.) Fig. 3 is a typical example, which 
shows how the magnitude of Bx- depends critically on the position of 
perigee, although the variation of Bx 2 , which also depends on other 
factors as previously stated, does not necessarily follow the same pat- 
tern as in Fig. 3 when the perigee is in other positions. At any rate, 
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Fig. 3 — B ± 2 vs longitude of perigee at start of orbit. 

the change of position of the perigee is a determining factor for the vari- 
ation of the day-average B~} and of the observed* spin-rate curve, as 
plotted in Fig. 2. The initial perigee of the Telstar satellite orbit (593 
miles in altitude) was north of the geomagnetic equator on the day of 
launch, as indicated in Fig. 1. The perigee advances in the direction 
of the orbital motion at a rate of approximately 2° per day. When the 
perigee was crossing the geomagnetic equator southward in the orbital 
plane, 7?T 2 first decreased and then increased, as indicated in the initial 
part of the TTj curve in Fig. 2. In the ascending part of the curve up 
to September 7, 1962, t = 60 days, more spin decay occiirred than would 
result from an exponential decay produced by a constant B ± 2 when the 
perigee was over the equator. This is why the corresponding part of 
the spin-rate plot is nearly a straight line instead of an exponential 
decay curve. From I = 60 days to t = 100 days, while the perigee was 
advancing northward toward the geomagnetic equator, BJ was level- 
ing off and then declining, resulting in an exponential decay as shown 
in the part of the spin-rate plot deviating from the extension of the 
straight line. From t = 100 days to t = 140 days, the perigee was en- 
tering the northern hemisphere, again getting into a stronger geomag- 
netic field indicated by the increasing BJ. As a result, this part of the 
spin-rate curve becomes nearly a straight line again, though of a dif- 

* The spin rate of the Telstar satellite was measured by J. S. Courtney-Pratt 
and his coworkers by means of the glint method (see Ref . 3 for details). It was 
also determined by C. C. Cutler and W. C. Jakes by way of measuring the fre- 
quency of the ripple in the amplitude of the radio signal received from the satel- 
lite. 
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ferent slope than the first one. Then the spin decay became exponen- 
tial again when the perigee was moving southward toward the equator 
from t = 140 to t = 190 days. All these indicate that the actual spin- 
rate curve would wiggle about a mean exponential curve as shown in 
Fig. 2. Note from Fig. 2 that BJ is lower when the perigee is in the 
southern hemisphere. This is due to the zone of depressed geomagnetic 
field strength previously mentioned. It is believed that, as the perigee 
keeps advancing in the orbital direction, B ± 2 (l) will continue to vary 
sinusoidally with time. Nevertheless, complete values of B ± 2 cannot be 
predicted for the entire useful life of the satellite because the spin axis 
changes its direction continuously because of perturbation of the elec- 
tromagnetic torques as well as occasional operation of the torque coil,* 
and because variation of the orbital parameters cannot be predicted 
accurately. Therefore, exact evaluation of the 1/e characteristic time of 
the nearly exponential spin decay cannot be obtained. Nevertheless, it 
may be approximately evaluated as follows. 

First, let us determine the value of p of the Telstar satellite from the 
observed spin-decay rate in the first 35 days, which is practically linear 
with time, withwo = 18.67 rad/sec (178.33 rpm)_aM = Oandw = 16.9 
rad/sec (161.2 rpm) at t = 35 days. Let us take B ± 2 and W as constants 
in (13) and, as t/r is small, we may expand exp (—t/r) up to the first- 
order term in t/r 



( l -¥')-£ 



u..,i i - '~r t) - :~y l 



Substituting into the above with w = 16.9, w = 18.67 rad/sec, I z = 
5.61 kg-m 2 , W = 1.6 X 10~ 7 joules, and t = 35 days = 3.024 X 10 fl 
sec, we find that 

pBJ = 1.748 X 10" 7 weberVohm. 

If we take B7 to be 1.572 X HT 10 weber 2 /meter 4 , then p = 1110 meter 4 / 
ohm. Now, we note from Fig. 2 that BJ varies between 0.0154 and 0.0190 
gauss 2 . If we take the average value of BJ for the entire useful life of the 
Telstar satellite, denoted as B ± 2 , to be 0.0177 gauss 2 ± 2 per cent, then 
the average 1/e characteristic time of the exponential spin decay with- 
out considering hysteresis losses is found to be 

r = -=L= = 330 days ± 2 per cent (15) 

pBJ 



* The torque coil consists of 200 turns of 32-gauge copper wire wound around the 
equator of the satellite. When current is turned on at a desired time, the magnetic 
moment of the coil will interact with the geomagnetic induction to produce torque 
for correction of the spin-axis direction. 
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for p = 1 110 meter /ohm. An exponential curve based on this mean 1/e 
time is plotted in Fig. 2; in addition, the actual spin-rate curve is also 
drawn (not to scale) in order to show that the latter curve is fluctuating 
about the mean exponential curve at a period of about 180 days. If the 
hysteresis losses are taken into account, the 1/e characteristic time is 
expressed as 

T h = t In [(«<, + W/2TpWT 2 )/(ao/e + W/2irpWT 2 )] (16) 

where W is the average value of W(t) for the entire useful life of the 
satellite. Let W be 1.6 ergs per cycle of rotation for a conservative esti- 
mate, then Th = 327 days ± 2 per cent. 

Based on the above range of the exponential decay rate, the satellite 
will spin at about 20 rpm at the end of two years from the day of launch. 
If the equatorial antennas had not been insulated, a separate calculation 
indicates that the spin rate after two years would be only about 3 rpm, 
which seems too low to insure attitude stabilization. 



III. SPIN PRECESSION AND PRECESSION DAMPING 

Before analyzing precession damping, let us first consider precessional 
motion of a spinning satellite produced by torques acting transversely 
to the spin axis. Suppose that the satellite, assumed here to be a rigid 
body, is initially spinning about its z-axis so that its initial angular 
momentum is J = 7 z o> = I z uz. When a transverse torque T is acting on 
the satellite for a time interval At, J changes to J by an amount A J, 
which is equal to the impulse TM, as shown in Fig. 4(a). The satellite 
will then perform precession with the spin axis, z, and the angular 
velocity, <o, no longer aligned with J. 

The torques causing precession consist of (a) gravitational torque, 
(b) eddy-current torque, and (c) torque of interaction between the 
residual magnetic dipole moment, M, and the geomagnetic field, H. 
The components of the gravitational torque, 3(gR Q 2 /p 3 )p X O-p (p = 
geocentric distance, g = gravitational acceleration at earth's surface, 
and R Q = earth's radius) normal to the spin axis are proportional to 
(Tz — Ix) or (I z — I u ) and are found to have a maximum value of 
0.65 X 10~ 6 ft-lb at perigee of the Telstar satellite orbit (based on the 
measured values of I x = .3.7140, /„ = 3.8252, and h = 4.1412 slug-ft 2 ). 
The eddy-current torque given in (6) can reach a maximum value of 
0.56 X 10" 6 ft-lb (for p, = 684 meter 4 /ohm and | B ]{ B ± | = 0.6 X 10" 10 
weber /meter 4 ). The maximum magnitude of the magnetic torque, 
M X H, is as high as 13.2 X 10" 6 ft-lb (for H = 31.84 amp-turns/m or 
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Fig. 4 — (a) Motion of Poinsot's inertia ellipsoid; (b) motion of oblate spher- 
oidal rigid body; and (c) motion of ball in the precession damper. 

2180 



SPIN DECAY 2181 

0.4 oersted and M = 0.502 X 10" 6 weber-meter*). Thus, both the 
gravitational and eddy-current torques are at least one order of magni- 
tude smaller than the magnetic dipole moment torque. 

For the analysis of precession damping, let us assume that J is tem- 
porarily an invariant, since T is so small (about 10~ 5 ft-lb maximum, as 
given above) that it takes a minimum time of about 1.5 days for J to 
change its direction by one degree, whereas the 1/e characteristic pre- 
cession damping time, t p , is only of the order of 30 minutes, as will be 
shown later. Thus, within the time interval comparable to t p the pre- 
cessional motion can be treated as torque-free. Such a motion can be 
pictured by Poinsot's geometrical construction (see Ref. 4, p. 161) in 
which the satellite's inertia ellipsoid rolls without slipping on the in- 
variant plane, which is a plane normal to J and tangent to the ellipsoid 
at a fixed distance from the origin of the ellipsoid [see Fig. 4(a)]. The 
curve traced out by the point of contact on the ellipsoid, known as the 
polhode, is the locus of the tip of <o in the body, while the curve on the 
invariant plane, known as the herpolhode, is the locus of the tip of <•> 
in the inertial space. To simplify the analysis of the precession damping, 
we further assume that the satellite is symmetric about its spin axis, or 
the inertia ellipsoid is an oblate spheroid with a transverse moment of 
inertia I(=I X = I v ) < h . In this case, the precession motion can be 
visualized, as shown in Fig. 4(b), as a body cone, zOw, rotating at an 
angular velocity, ii, on an immovable space cone, JOu, which is rotating 
at an angular velocity, o>i , along the fixed direction of J. The line of 
tangency between these two cones is the instantaneous axis of rotation 
of the body or the angular velocity, to, which is the sum of £1 and on . 
The analytic solution of such a torque-free precession motion of an 
oblate spheroidal body is obtained (see Ref. 4, p. 162) for the angular 
velocity to = u x x + w u y + u 2 z expressed in the body coordinates with 
components 

Wj = a>x sin tit 

u u = Wj. cos tot (17) 

o>z = w.( = constant) 



where 



C-9 4 ) 



n = [—t^ «« (is) 



* Meusured by M. S. Glass and 1). P. Brady. 
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and, as is apparent from Fig. 4(b) 

Wi = (w z — ft) tan 6 — -f- o) z tan 6 = constant. (19) 

Here, it is obvious that the precession angle, 0, between the spin axis, z, 
and J is a constant, providing there exists no precession energy dissipa- 
tion. 

The precession energy is the difference between two energy states with 
and without precession. We use the same assumption as in Section I: 
that the precession energy dissipation produces negligibly small torque 
to the rigid body motion. Thus, during the time interval when the 
precession is substantially damped out, the angular momentum can be 
treated as an invariant. The kinetic energy in the presence of precession 
is 

E = \Joix + \I z Uz 

where co x 2 = (*h? + <*>/, and as shown in Section I the minimum energy 
state occurs when the precession is completely damped out, i.e. 

where co/ can be found from the invariant angular momentum 

,2 T 2 2 . r 2 2 r 2 II 

J = iWi ■+■ lz 0i z = 1 z l»z 

or 

Therefore, the precession energy is 



E p = E - E min = - 



K-i)'-- 

or, by virtue of (19) 

E p = i (j - 1) W tan 2 6. 

Note that when / = 7, for a spherical satellite there is no precession 
energy. 

To dissipate the precession energy, the satellite is equipped with a 
pair of curved aluminum tubes filled with neon gas at one atmosphere, 



(20) 
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each containing a tungsten ball of radius r and mass m [see Fig. 4(c)]. 
The ball is slightly smaller than the inside diameter of the tube, and the 
curved tubes are installed concavely towards the spin axis with their 
bisecting radii of curvature, R, perpendicular to the spin axis at the 
center of mass of the satellite. When the satellite is rotating precisely 
about its spin axis, the balls are stationary at the middle of the tubes. 
However, when precession occurs — i.e., when w is processing along the 
body cone or following the polhode in the inertia ellipsoid — the balls are 
forced to move back and forth against the viscous friction of the gas as 
well as the inviscid friction between the balls and the tubes. Hence, the 
precession energy is dissipated into heat through the resistances to the 
motion of the balls, resulting in the attenuation of the precession angle, 
0, or in the realignment of the spin axis, z, and <o with J. 

To derive the equations of motion of the ball, we assume that the ball 
rolls on the tube without sliding. The equation of the ball's rotational 
motion about its center of mass can be immediately written in terms of 
the angle, a, from the ?/-axis (the tubes are assumed to lie in the yz- 
plane) 

?™' 2 (f«) =fr-N (21) 

where / is the force acting at the point of contact, and N the resistance 
moment due to rolling friction. The position vector of the center of mass 
of the ball, as shown in Fig. 4(c), is given in the body coordinates as 

D = [R(l - cos a) - b]y + R sin az. (22) 

The equation of the translational motion of the ball is then 

m -r#--q = — / — cRa (23) 

where q = sin ay + cos az is the tangential unit vector in the direction 
of a and c the coefficient of viscous friction. In performing the differentia- 
tion of D with respect to time, one should be aware of the fact that the 
angular velocity, to, as given in (17) in the body coordinates of a pro- 
cessing body, is changing in direction in an inertial space and its time 
derivative is to = to X ft, where ft = ttz. Therefore, it should be noted, 
for example, that (d/dt)y = to X y and (d 2 /dt 2 )y = (to X ft) X y + 
^ X (to X y). Upon differentiating D in (22) twice with respect to time, 
substituting into (23), and using (21) to eliminate /, we obtain a 
nonlinear second-order equation for a 
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a — — = o) z \R(l — cos a) — b] sin a + - a>„(to z — J2) sin 2 a 
lit 7 

5 5 

— =o» ± sin a cos a + — w w (w g + fl)[fl(l — cos a) — b] cos a (24) 

5c . a 5N 



a — 



lm | a | 7mrR 

where the sign of N has been chosen in such a way as to make the resist- 
ance moment always oppose the rotational motion of the ball, and o) y = 
«! cos Q£ = (o) z — £2) tan 9 cos 9,t. As the radius of curvature, R( = 15 ft), 
of the tube used on the Telstar satellite is much larger than its length, 
L(ph1A ft), the maximum angle of a is very small, viz., a m = L/2R = 
0.0465 radian «1, where a m is the subtended half angle of the curved 
tube. In order for the balls to move back and forth without bottoming 
with the ends of the tubes, the precession angle, 6, should be of the same 
order as a. Thus, for small a and 6, (24) can be linearized to the following 

a + 2na + P 2 a + j^-. K = q cos fii (25) 

I <* I 

where 2m = 5c/7m, P 2 = 5bu,*/7R, K = 5N/7mrR, and q = (5b/7R) • 
6{<j) z — ii 2 ). Because of the presence of the rolling friction term, the 
above equation can be solved only for each half cycle. 

An experiment has been conducted by the author to determine the 
resistance moment, N, that the tungsten ball encounters when rolling 
on the aluminum tube. For convenience, N is expressed in terms of a 
resistance force, F, acting at the center of the ball: i.e., N = Fr, and F is 
determined experimentally to be F = 0.0002 lb. In another experiment 
devised to measure the coefficient of viscous friction, c, for the 0.484- 
inch ball moving in the tube (nominal inside diameter = 0.495 in.) filled 
with neon gas at one atmosphere, it is found that c = 0.00193 lb-sec/ft. 
The ratio of the energy dissipation per cycle due to the viscous friction 
(in the steady-state forced oscillation case) and that due to the rolling 
friction can be shown to be E,/E, = | wcRQa m /4F |. With a m = 0.0465 
radian, R = 15 ft, 1 12 [ = 1 .50 rad/sec, corresponding to the case of the 
Telstar satellite at 178.33 rpm, the above ratio is about 8. This indi- 
cates that energy dissipation per cycle due to rolling friction is approxi- 
mately one order of magnitude smaller than that due to viscous friction 
at the indicated spin rate. If the rolling friction term is neglected, (25) 
becomes 

a + 2na + P 2 a = q cos Qt (26) 
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for which the steady-state forced oscillation solution is 

a = a COS (fl« - j8) (27) 

where 



«„ = { 1 - - 
Oil 



P 2 , 



'V 4nVT 
V + P 4 J 

cos jS = (P 2 - 12 2 )[(P 2 - ft 2 ) 2 + 4nV]- J (28) 

sin/3 = 2J2n[(P 2 - ft 2 ) 2 + 4n 2 fi 2 ]~\ 

The energy dissipation due to viscous friction per period of oscillation of 
the ball (or per period of precession) is 

E v = / ctiWdt = cRWSl / sin 2 (Sit - 0) d(Qt) = vcRWQ. (29) 

•'0 "O 

The time-average rate of energy dissipation per period of precession 
appears to be 

dE v E v 1 p2 2 2 



I) 2 



Equating the negative of the above to the rate of change of the preces- 
sional energy, E p , in (20) for the case of a small angle, tan 6 « 0, an 
equation for the change of the precession angle is obtained 

L A -r 2„, 1 r.2 rt 2 



j - i J i^-ee = -JcTO. (31) 

Substituting into the above with ur given by (28) and integrating, we 
obtain an exponential decay of with time 

6 = o e~' /Tp (32) 

where d = d(t = 0) and r p i.s the characteristic time given as 

T " = 7nm#(X -Y)X*(2 - X) 2 LV " ?7 + ~^~J (33) 

with X = 7-/7 ( > 1 ) . If the time-average rate of energy dissipation per 
cycle due to rolling friction, ±FRa 1 ft |/2tt, is included in (31), the 
solution for becomes 

6 = (0 O + k)e~ t,r ' - h (34) 
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where 



10F 

h = 



^[(-^v-fl- ») 



7irmnRio z (\ - 1)X(2 
The l/e characteristic time in this case is 

'i + f 

r.' = r. In | 5. | (36) 

which involves the initial angle, 6 . Because of the rolling friction, the 
precession will be damped out within a finite time, i.e., 8 = at t = 
t p In (1 + $o/h). A numerical computation, corresponding to the 
Telstar satellite physical constants, shows that h is considerably smaller 
than one degree. This indicates that if the precession angle, 0, is substan- 
tially greater than one degree, the l/e characteristic time given in (33) 
for viscous friction alone should be used for convenience, since it is 
independent of the initial condition. 

From (26) it is seen that if the natural frequency, P, is made equal to 
the frequency of the forcing term, ft, i.e., if 

R = w^w b (37) 

then the oscillating motion of the ball is hi resonance with the precession 
motion of the satellite. As a consequence, the l/e characteristic time 
becomes much shorter 

28n/ z (X - 1) , . 

Tpr 5mPXW(2 -X) 2 ' K } 

Unfortunately, such a tuned damper cannot be obtained for the Telstar 
satellite, since the ratio of moments of inertia (I/I z = 0.897, 0.925 or 
X = 1.114, 1.08) is close to unity, and as the tubes are placed outside 
of the electronics package, 6 cannot be made too small. Therefore, in 
view of (37), a tuned damper for the Telstar satellite would have to be 
of an exceedingly large radius of curvature (R = 57-117 ft for b = 
1.046 ft). In this case, the tubes become practically straight, and the 
motion of the balls may not necessarily be at resonance with the pre- 
cession of the satellite. The equation of motion of a ball in a straight tube 
can be easily obtained by multiplying (26) through with R and then 
letting R — > oo or P — > 0; in a similar way, the l/e characteristic time 
can be obtained. Nevertheless, such a straight or nearly straight tube 
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damper will not be used for the essential reason that, in case of its 
misalignment with the spin axis, no damping whatsoever will be obtained 
when the precession angle is smaller than the misalignment angle. Another 
type of damper can be obtained if the curved tube is placed concavely 
away from the spin axis. The equation of motion of the ball in such a tube 
can be easily shown to be the same as (26) except that the sign of the a 
term is negative, hence forming an unstable system. The ball, which 
rests at one end of the tube, will move toward the other end at a high 
speed when a component of to is tangent to the former end of the tube 
and is large enough to make the centrifugal force overcome the static 
friction. The ball will move back and forth twice in each period of pre- 
cession. Let the kinetic energy of the ball when it reaches the other end 
of the tube equal the centrifugal force times the distance traveled by the 
ball perpendicular to to. If we assume that the kinetic energy of the ball 
is completely absorbed by bottoming at each end of the tube, it can be 
shown that the decay of the precession angle is parabolic with time. 
Such a damper can effectively reduce the precession even when the ratio 
I/I z is close to unity, but it will not damp out a precession angle less than 
about three degrees; thus, it was not adopted for use with the Telstar 
satellite. 

After comparing the advantages and disadvantages of the several 
dampers discussed, the untuned concave damper shown in Figs. 4(b) and 
(c) was finally chosen for the Telstar satellite, although its damping time 
is somewhat larger than that of the others. This choice was made because 
the theoretical \/c characteristic damping time, given in (33), is calcu- 
lated to be a maximum of about three minutes for a ratio of I/I 2 up to 
0.95, a spin-rate range of 20-180 rpm, and for the parameters given be- 
low. Such a damping time is acceptable even if it is one order of magni- 
tude larger, in view of the slow rate of change of the angular momentum 
due to the small transverse torques previously calculated. The chosen 
parameters are: R = 15 ft (a large value, although the tube still has 
noticeable curvature), m = 0.0021 slug (for two tungsten balls of 0.484 
in. diameter; tungsten is chosen for its high density), and n = 0.65 sec -1 
or c = 0.00103 lb-sec/ft (corresponding to neon gas, which is chosen for 
its high viscosity). [Note: For a tuned damper with 1/e damping time 
as given in (38), a gas with low viscosity should preferably be used.] 

It is necessary that the theoretical 1/e characteristic time t p of the 
untuned concave damper should be compared with experimental results 
for the following reasons. Formula (33) is obtained from the linearized 
analysis of the motion of the ball under the assumptions of small ampli- 
tude of the motion and an axisymmetric spinning body. In the actual 
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case, the Telstar satellite is dynamically not axisymmetric. Also, as the 
tube is limited in length, bottoming will occur when the precession 
angle is larger than about 3.5°; the motion of the ball will then be dis- 
turbed and will not follow (26). An experiment has been devised by G. T. 
Kossyk which consists of an air-bearing supported flywheel (I z = 5.14, 
/„.ux = 4.675, I min = 4.404 slug-ft 2 ) mounted with two precession dam- 
pers as shown in Fig. 5. The flywheel is driven to reach a certain initial 
speed about a skew axis making a desired angle with the axis of sym- 
metry, which is the principal axis of maximum moment of inertia. As 
soon as the drive is released, the spinning flywheel performs a preces- 



PRECESSION 
DAMPER 




Fig. 5 — Schematic layout of precession damping experiment. 
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sional motion with known initial angular speed and initial precession 
angle. The decay of the precession angle is recorded through an optical 
tracking device for two different cases, with and without the balls in 
the damper tubes. The difference between these two recordings is the net 
effect due to the precession damper, excluding all other effects due to air 
resistance, gravity, etc. The balls were observed to be moving, and bot- 
toming was clearly heard. The experimentally determined r p is found to 
be about four times larger than the theoretical t p calculated on 
the flywheel based on the mean value of the transverse moments of in- 
ertia and about nine times larger based on the minimum transverse 
moment of inertia. Although the Telstar satellite has different moments 
of inertia from those of the flywheel and a higher I max /I z ratio, it is 
believed that the actual r p should not be greater than the theoretical t p 
in (33) by more than one order of magnitude. 

For a conservative estimate of the precession damping time of the 
Telstar satellite, let us multiply (33) by a factor of nine and use the 
following physical constants: h = 4.1412, I max = 3.8252 slug-ft 2 , X min = 
/.//max = 1-08, m = 0.0021 slug, n = 0.65 sec" 1 , R = 15 ft, 6 = 1.046 
ft. Equation (33) is then reduced to 

t„ = 19 (0.76 + ^ ) minutes 



'i' 



H+f: 



which is relatively independent of the spin rate in the range of interest. 
At 178.33, 65, and 24 rpm, the maximum 1/e characteristic precession 
damping times are 14.8, 16.2, and 27.6 minutes, respectively. 

IV. DRIFT OF SPIN AXIS 

We have shown in the Section III that the major transverse torque 
causing spin precession is contributed by the residual magnetic dipole 
moment along the spin axis, M. (M is found to be pointing toward the 
rocket-mount end of the Telstar satellite.) The torque produced by the 
residual magnetic dipole moment normal to the spin axis is mostly 
averaged out because of the spinning motion ; other transverse torques, 
produced by eddy currents and gravity, are all one order of magnitude 
smaller than that produced by M, as shown previously. Therefore, in the 
qualitative analysis of the spin-axis drift in this section it is sufficient to 
take only M into account. 

The initial direction of the Telstar satellite spin axis on the day of 
launch was 82.3° right ascension and —65.6° declination, as represented 
by the initial angular momentum, J , (see Fig. 1) in the nonrotating 
coordinate system O-XYZ, where OX points toward the vernal equinox 
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and OZ is along the earth's spin axis. The geomagnetic field, as shown 
schematically in Fig. 1, is rotating about OZ at the earth's spin rate. 
When the satellite is traveling along its orbit, it finds that the geomag- 
netic induction, B, generates a nearly conical surface with respect to 
Jo or to O-XYZ, with the axis of the cone pointing in the direction of the 
orbital angular momentum J orb (see Fig. 6). Let us pass a plane through 
O normal to J , project the conical surface (or B) onto the plane, and 
construct the time-average resultant of the projection, Bj. . Then B x 
will interact with the magnetic dipole moment M (pointing in the nega- 
tive direction of J or w) to produce a transverse torque, T x , which 
causes spin precession. Because of the rotation of the geomagnetic field, 
which has anomalies, and because of the apsidal advance in the orbit, 
the conical surface generated by B has a different area every orbit. 







Fig. 6 — Schematic illustration of spin-axis drift. 
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Furthermore, due to the precession or the nodal regression of the orbital 
plane, the cone gradually changes its direction in the O-XYZ coordinates, 
indicated in Fig. 6 as the rotation of J or b about OZ. As a result, Bj. 
and hence T A gradually change in both direction and magnitude, while 
the angular momentum, J, of the satellite follows the pattern of varia- 
tion of Tj. , as shown in Fig. G. Because of spin precession, the spin axis 
is turning around J, yet does not deviate away from J as a consequence 
of precession damping. Therefore, the tip of the spin axis describes a 
spiral path, shown in an exaggerated way in Fig. 6. Such a phenomenon 
is termed the drift of the spin axis. The pattern* of the drift is deter- 
mined by the orbit and by the orientation of M, whereas the rate of 
drift depends on the orbit and the magnitude of M. 

A rough estimate of the rate of drift of the Telstar satellite spin axis 
can be given. Let us assume that the transverse torque, Tj. , keeps acting 
on the satellite perpendicular to the angular momentum, J, despite the 
fact that J continuously changes its direction as time goes on. This 
assumption is justified by the fact that the precession dampers work 
properly, so that the spin axis is virtually in line with J. Let us disregard 
the retarding torque at this point. Then, from the principle of angular 
momentum about the center of mass of the satellite 

2 = T A (39) 

at 

we find that after a time interval At the angular momentum changes to a 
new position by an angle 

M = T -^. (40) 

In the above we have kept J = I z <a, or o> at a constant magnitude, be- 
cause we have not considered the retarding torque. To evaluate Ad in a 
time interval of one week, let us substitute into the above with At = 
6.048 X 10 5 sec, J, = 5.61 kg-m 2 , T x = MH ± , where M = 0.562 X 
10~ 6 weber-meter and H ± = H cos y (H = time-average value of the 
geomagnetic field on the Telstar satellite orbit and H x is the com- 
ponent of H along the spin axis; y is the angle between H and the spin 
axis). Then (40) is reduced numerically to 

n -.H cos 7 , 
Ad = 274 - degree 

CO 

where H is in oersteds and w in rad/sec. If H is taken to be 0.2 oersted 



* For details of the pattern and rate of drift, see Ref. 5. 
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and 7 to be 45°, then at the initial spin rate of 178.33 rpm or u = 18.67 
rad/sec, we find A0 = 2.1° per week, which is very close to the observed 
value (approximately 2° per week). This drift rate should increase 
exponentially with time because of the exponential decay of oj resulting 
from the action of the retarding torque. It appears from the above esti- 
mate that the satellite's residual magnetic dipole moment along the spin 
axis did not change drastically due to launching. 

V. SUMMARY AND CONCLUSIONS 

The dynamics problems for the spin-stabilized Telstar satellite, char- 
acterized by spin decay, spin-precession damping, and spin-axis drift, 
have been studied in this paper. In the section on spin decay, the nature 
of the retarding torque due to eddy-current losses has been analyzed. 
The observed decay phenomena are largely explained from the computed 
B± 2 , taking into account the anomalies of the geomagnetic field, the 
variations of orbital parameters, and the change of the spin-axis direc- 
tion. The 1/e characteristic time of the nearly exponential spin decay is 
estimated to be about 330 days ±2 per cent by extrapolation from the 
observed data. This indicates that at the end of two years from the day 
of launch the Telstar satellite will spin at approximately 20 rpm. It is 
believed that motion at such a spin rate is still relatively stable with 
respect to precession. 

For the spin precession, it is found that the transverse torque is pro- 
duced mainly by the residual magnetic dipole moment along the spin 
axis. The precessional motion of a spinning satellite is illustrated by 
means of Poinsot's geometrical constructions. A few different types of 
precession dampers have been considered. Linear analysis of the motion 
of the ball in the concave type damper has been made, from which ex- 
plicit expression of the theoretical 1/e characteristic precession damping 
time is obtained. Based on the analysis, it was possible to make a proper 
design of the damper. An experimental comparison of the theoretical 
1/e time enables us to estimate the actual 1/e time to be about 30 min- 
utes maximum. It is concluded that this damping time is acceptable for 
the computed magnitude of the transverse torques, and in fact, no pre- 
cession angle larger than 0.5° has yet been observed on the Telstar 
satellite. 

In discussing the problem of spin-axis drift, only a brief qualitative 
description is given to illustrate the fundamental mechanism; also, an 
approximate quantitative analysis is shown for an order-of-magnitude 
estimate of the drift rate. The observed continuous drift of the spin axis 
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of the Telstar satellite is evidence of proper functioning of the precession 
dampers. 

The above three problems, which are caused essentially by electro- 
magnetic torques, can be summarized into one of the important dynamics 
design criteria of a spin-stabilized satellite: i.e., evaluation of the maxi- 
mum allowable eddy-current losses and residual magnetic dipole moment 
for specified useful life and orbit of the satellite. Other basic dynamical 
requirements are worth remarking here. The spin axis should necessarily 
have a maximum moment of inertia because of provision of precessional 
energy dissipation and because of elastic energy dissipation, since the 
satellite is not a perfectly rigid body. This moment of inertia should also 
be made as large as possible for a given weight and size of the satellite, 
in order to make the satellite more stable and to increase the lifetime for 
the same initial spin rate. Furthermore, the ratio of the moment of in- 
ertia about the spin axis to those about the transverse axes should be 
made large enough to yield an adequate precession damping time. 
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